Some exact eigenstates of the generalized squeezing operator 
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We explicitly construct the states that are invariant under the action of the generalized squeezing 
operator exp (za^ — z*a!') with integer fc > 3. There are k such states for a given value of k. We 
show that in the number representation the states behave as n"'"''* for large n's. This implies that 
that for any > 3 the states are normalizable. However, for a given k the expectation values of 
operators of the form (a^a)-' diverge for the integer j > (fc/2 — 1). For fc = 3 we also give an explicit 
form of these states in the momentum representation in terms of Bessel functions. 



I. INTRODUCTION 

The concept of squeezing plays one of the central roles 
in quantum optics. Review articles [TJ [5] provide an ex- 
tensive hst of references that deal with various aspects 
of squeezing. The mathematical realization of squeezed 
state in the simplest case is given in terms of the squeez- 
ing operator S{z) = exp((za^^ — z*a?)/2) acting on the 
vacuum state. Over time, attempts to generalize this op- 
erator to include higher order processes have been made. 
It turned out, however, that the so-called naive general- 
izations of the squeezing operator of the form 



Uk{z) — exp (za^ 



(1) 



with integer A; > 3 proved to be quite nontrivial. On one 
hand, it was shown by Fisher et al. [3] that the vacuum 
to vacuum probability amplitude (0|/7fc(z)|0) has a zero 
radius of convergence as a power series with respect to z, 
for k > 2. On the other hand, Braunstein and McLach- 
lan demonstrated numerically [H [5] that such expressions 
can still be well defined. Some properties of operators 
Uk{z) were discussed by Nagel [6 . In Ref.[7 Elyutin and 
Klyshko showed that the average occupation number of 
an arbitrary initial state acted upon by the unitary oper- 
ator Uji^z) diverges to infinity for a finite z. All of these 
challenges lead to various attempts to modify the opera- 
tors and a in Uk{z) in such a way that no question of 
convergence would arise [8l[9l[10]. In this paper we con- 
struct k exact degenerate eigenstates of operator Uk{z) 
that have an eigenvalue equal to one and investigate some 
of their properties. 



II. INVARIANT STATES OF OPERATOR Uk(z) 
IN THE NUMBER REPRESENTATION 



Consider operator za^ — z*a*' that appears in the ex- 
ponent of Uk{z). Its eigenvalues are not known except for 
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the cases of fc = 1,2 [6,. In general, the eigenvalue prob- 
lem of this operator written in the number representation 
leads to a three term recurrence relation. However, in the 
case of zero eigenvalue (assuming it exists) the recurrence 
relation becomes a two term one and, because of this 
simplification, explicit eigenstates can be obtained. It is 



clear that since these eigenstates of za 



have zero 



eigenvalues they are also eigenstates of Uk{z) with eigen- 
value one, or , in other words, invariant under the action 
of Uk{z). Without the loss of generality we can limit 
ourselves to the case where z is imaginary z = — = ir. 

Thus, we consider zero eigenstates of the Hermitian 
operator 



(2) 



Even though we will be primarily interested in the case 
of fc > 3, results of this section also apply to fc = 1,2. 
Note that operator Uk{z) commutes with operator Gk ~ 
exjp{i^a'^a). 

In the case of zero eigenvalue the eigenvalue equation 
in the number representation reads 



{n\H,\i,^) ^ 0. 



(3) 



where a is the degeneracy index. Acting on {n\ from the 
left with the creation and annihilation operators in Hk 
we obtain the following recurrence relation 



(n + fc|Vn 



ln{n- l){n- 2)...(n- fc+ 1) 
{n + l){n + 2)...{n + k) 



n-k\rl^^) 



(4) 

as well as the following fc conditions 

(fclV-fe") - 0, (fc + m) = 0, ...(2fc - llV-n = 0. (5) 

By iteratively applying recurrence relation of Eq. ^ 
starting with (0|V'fc), {l\ip'^), ...(fc — 11?/;^), we obtain after 
some algebra fc degenerate zero eigenstates of Hk- The 
non-vanishing components of these eigenstates have the 
following form 



(a-h2mfc|V'fe) 



(2fc) 



km 



\/c{k,a) ^(a + 2mfc)! 



a + i 



2k 



bittner6un.edu 



(6) 



2 



Here a is the degeneracy index that runs from to (fc— 1) 
and integer to (not to be confused with the occupation 
number) runs from to infinity. V(x) denotes the gamma 



function and c(fc,Q!) is the normalization constant given 

by 



c(fc,a) =fe+iFfc 1 



1 



2fc ' 2k 



2k 



2k 



k + 2- 



2fc + a, \ n- (r(^))' 



2k 



2k 



(7) 
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FIG. 1: Functions {nlip'^) (triangles), {n\jpl) (stars), and 
{n\ip3) (diamonds). Lines are added to lead the eye. 



where pFq{xi,X2, ...,Xp;yi,y2, ■■■,yq; z) is the generaUzed 
hypergeometric series. Eigenfunctions of Eq. ^ are 
monotonically decreasing functions of m. Fig. 1 shows 
the graphs for {n\ip^)'s for k — 3. Note each of |'0fc)'s 
is also an eigenstates of operator Gk with eigenvalue 
exp(i^Q;). The asymptotic behavior of functions (|6| 
for large m can be obtained with the help of Stirling s 
expansions for the gamma functions and the factorial in 
Eq. ([6|. After some algebra we obtain for large to's 



(a + 2mfc|V') d{k, a)m~^/'^. 
Here the pref actor d{k,a) is given by 

1 (27r)(^^-i)/^ 



d{k, a) 



(8) 



(9) 



Since the square of the eigenstate behaves as mT^^'^ for 
large to we can conclude that the norm is finite for any 
fc > 3. For A; = 1 and 2 the norm diverges in agreement 



with the known exact results for these cases [6] . Similarly, 
we can see that the average for the number operator h = 
a^a is divergent for k < 5, diverges for k < 7, etc. In 
general, for a given k the expectation values of operators 
of the form diverge for the integer j>(fc/2 — 1). 

If we define dimensionless coordinate and momentum 
operators x = \J\{^'^ + a)i P = i\J\{0'^ ~ then their 
expectation values for states in Eq. ([6| vanish. However, 
if superpositions of the degenerate states ^ are consid- 
ered then, in general, the average of x and p will diverge 
for fc < 3 but converge for k > A. Expectation values 
of x"^ and behave in the same way as that for a^a, 
namely, diverge for fc = 3, 4, but remain finite for fc > 5. 

Divergence of the expectation value of the number op- 
erator for fc = 3, 4 implies infinite average energy for these 
states. Thus, they cannot be produced experimentally. 
The situation is more subtle for fc > 5 since in this case 
the average energy is finite. 



III. INVARIANT STATES OF U^iz) IN THE 
MOMENTUM REPRESENTATION 



Since eigenstates decay slowly as functions of 

n, it is of interest to consider their behavior in a contin- 
uum basis, such as coordinate or momentum representa- 
tions. In this section we will consider a specific case of 
momentum representation for fc = 3. We chose momen- 
tum over coordinate representation to demonstrate an 
interesting mathematical point that will be mentioned 
below. Rewriting a and a) in terms of dimensionless co- 
ordinate and momentum operators as a = ^J~^{x+^p) and 

= \J\{x — ip) and inserting them into we obtain 
the following eigenvalue equation for the zero eigenvalue 



dp^ 



dp 



1 



(10) 



Here we suppress subscript 3 for fc = 3 to simplify the 
notation. This is a second order differential equation and 
its two independent solutions are given by jllj 



- (275) ■ 
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FIG. 2: Functions (p|<^^) = 
{p|^2) = iVo(^) (dashes). 



Jo (^) (solid line) 



and 



2V3 



(11) 



where 0/0(2;) ^^id iVo(x) are the zeroth order Bessel func- 
tions of the first and second kind, respectively. Note that 
functions {p\ip^) and (p|</3^) are neither orthogonal to each 
other nor normalized. Plots of these functions are shown 
on Fig. [2] The behavior of {p\(p^) and (p|</?^) for large p 
is determined by the asymptotic behavior of the Bessel 
functions 



2(31/4) / p2 ^ 

— = — cos — = - - 
V^P V2\/3 4 

2(31/4) . / p2 ^ 

— = — sm — - - - 
\^P \2V3 4 



(12) 
(13) 



Function {p\ip'^) has a logarithmic singularity at p = 0. 
Both {p\(p^) and {p\(p'^) are even functions of p and, there- 



fore must be linear combinations of eigenstates 
and {p\ip'^). The obvious question then is what hap- 
pened to the third function \ip^) which must be odd in 
the momentum representation. The answer to this ques- 
tion comes from noting that Eq. 10 is singular at p = 0. 

Toll 



This becomes obvious once both sides of Eq. ( |iup are 
divided over by p to bring the equation to the standard 
form. The third solution is obtained by reflecting {p\ip^) 
taken from —00 to with respect to the p axis. Thus, 



Jq 



(P), 



(14) 



where sgn(p) is the sign function. It is e asy to verify that 
{p\(p^) is indeed a solution of Eq. (10 1 since the differ- 



entiations of at the vicinity of the "step" at give zero 
contribution. Thus, for the case of fc = 3 eigenfunctions 
(p|V'")'s in Eq. ^ are expressed through (pl^p^), {p\^'^), 
and (pI'^^) as follows 



(p|V°) = ao(p|/)+5o(p|^2) 



(15) 
(16) 
(17) 



Coefficients oq , bo, 02 , and 62 can be found and expressed 
through known functions in various ways. In particular, 
they can be expressed through the overlaps of\(p^), |</3^), 
Itp^), and with the ground and second excited states 
of the harmonic oscillator as follows 



ao 

0.2 



D 

D 



bo^mp^^ (18) 



D 



where 



and 



D = (0|^i)(2|^2)-(0|^2)(2|vpi) 



, 2(31/4) /2-^/3^ 



(19) 



(20) 



(21) 



(31/4 ^ 33/4) ^ / 2- _ 272(3^/4) ^ / 2-^/3 



7r3/4 
1 



\/2(33/4)7r5/4 



r3/4 



(22) 
(23) 
(24) 
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(25) 
(26) 
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FIG. 3: Functions {p\ip°) (thick solid line), -i{p\ijj'^) (dashed 
line), and {p\ip^) (thin solid line). 



resentation - it is solvable in the Sturm-Liouville sense, 
namely, if we define functions 

/+(p) = Jo ('^) 9{p)e{i), 



V2V3, 

frip) = (2^) 



(28) 



where 0(x) is the theta function, it can be verified that 



(1-P) 



p'fHp)- (29) 



Thus, fj^{p) is an eigenfunction with eigenvalue |(1 — Z^) 
and weight function p^. Functions fj^{p) form a complete 
set, 



1 /"^ 1 

^ / dii{f+ip)f+{p') + frip)frip')) = -,5{p-p'). 

(30) 

It appears, however, that this solution cannot be used to 
the spectrum of the exponential operator U'i{z). 



K{x) and E{x) are the complete elliptic integrals of the 
first and second kind, respectively, and F{a, (3; 7; x) — 
2 Fi(a, /3; 7; x) is Gauss's hypergeometric function. Co- 
efficient fli in Eq. (16) is given by 



%1T 



ai 



l/4p( 



6r(^)r(i) 



(27) 



Functions (pIV'") are shown on Fig. Their asymp- 

totic behavior is essentially determined by the a symp- 
totics of {p\ip^) and {p\'p'^) given in Eqs. ( 12|l3 1. All 
(p|'0")'s show slow oscillating decay for large (pIV"") 
and (pIV'^) have logarithmic singularities at p = 0. All 
three functions, however, are square integrable in agree- 
ment with the results of Sec. |IT] Note that explicit solu- 
tions can also be obtained in the coordinate representa- 
tion either by Fourier transforming (p|'0")'s or by solving 
Eq. (10 1 rewritten in the coordinate representation. We 
will not consider them in this paper. Note, however, that 



in the coordinate representation Eq. ( 10 ) is the third or 



der differential equation and, therefore, the issue of the 
"missing solution" does not arise. 

Finally, let us note the following interesting property 
of the eigenvalue equation for in the momentum rep- 



IV. CONCLUSION 

We explicitly constructed some of the eigenstates 
of the generalized squeezing operator Uk{z) — 

exp (za^ — z*a'') in the number representations and 
showed that they are normalizable for fc > 3 but have 
divergent expectation values for operators (a^a)^ for the 
integer j > (fc/2 — 1). We produced only k eigenstates 

of Uk{z) ~ exp (za^^ — z*a'^). If we assume that the 
remaining eigenstates of Uk{z) have similar convergence 
properties this would imply that operator Uk{z) has a 
spectral resolution in the Hilbert space. Moreover, the 
states that we found can become useful for approximate 
treatments of operator Uk{z). In particular, if these op- 
erators are approximated by finite dimensional matrices, 
then the suitable basis sets can be chosen to have con- 
vergence properties similar to the states that were con- 
sidered in this paper. Due to their interesting properties 
the states can be of also interest for mathematical 
physics applications, i. e. in such areas as the theory of 
generalized squeezed and coherent states [T] or the theory 
of wavelets [12]. 
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